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abstract 


Rapid  bioassessment  data  were  used  to 
evaluate  Rosgen  stream  npe  designation  and 
several  other  stream  classification  vanables  as 
methods  for  partitioning  variation  among  the 
benthic  macromvertebrate  communities  of  high- 
quality,  cold  water  fisheries  at  the  stream  reach 
spatial  scale.  Benthic  macroinvertebrate 
communities  were  characterized  using  five 
metncs:  (1)  taxonomic  richness,  (2)  modified 

Ephemeroptera.  Placoptera,  Trichoptera  (EPT) 
Index,  (3)  modified  Hilsenhoff  Biotic  Index, 
(4)  percent  dominant  taxa,  and  (5)  modified 
percent  mayfly  index.  No  sigmficant  differences 
were  observed  in  the  metric  scores  of  sample  sites 
based  on  stream  order,  drainage  area  class,  surface 
water  slope  class,  or  Rosgen  Level  I stream  type 
designation.  Although  Rosgen  Level  II 
classification  resulted  in  excessive  fragmentation 
of  the  data  set,  significant  differences  w-ere 
observed  in  the  taxa  richness  and  modified 
percent  mayfly  index  scores  of  some  of  the 
Level  II  stream  tvpes.  How^ever.  these  differences 
were  more  readily  attributed  to  ecoregion 
designation  than  the  stream  reach  channel 
morphology  features  used  in  Level  II 
classification. 

A third  stream  classification  variable  based  on 
stream  reach  morphology  was  developed  dunng 
this  study.  This  classification  variable,  termed  the 
"modified  Rosgen"  stream  type,  involved 
classifv'ing  sample  sites  based  on  their  Rosgen 
Level  II  designation  without  the  surface  water 
slope  subcategor\-.  Significant  differences  were 
observed  in  the  modified  EPT  and  the  percent 
dominant  taxa  index  scores  of  some  of  the 
modified  Rosgen  stream  types.  The  modified 
Rosgen  stream  type  designation  was  the  most 


effective  stream  reach  classification  vanable  for 
partitioning  variation  among  the  benthic 
macromvertebrate  communities  of  the  high- 
quality',  cold  water  fishenes  (HQ-CWFs)  included 
in  this  study. 

INTRODUCTION 

In  an  effort  to  comply  with  federal 
antidegradation  requirements,  the  Pennsylvania 
Department  of  Enwonmental  Protection  (Pa. 
DEP)  developed  the  Pennsylvania  Special 
Protection  Waters  Program  (Shertzer,  1992).  The 
program,  which  includes  streams  or  watersheds 
that  exhibit  exceptional  water  quality  and  other 
environmental  features,  protects  selected 
waterways  by  establishing  higher  water  qualitv 
standards  than  those  that  typically  apply  to  the 
surface  waters  of  the  commonwealth.  Streams 
and  rivers  that  are  being  considered  for  inclusion 
in  the  program  usually  undergo  an  assessment  of 
their  "ecological  significance,"  which  is 
determined  based  on  their  EPA  Rapid 
Bioassessment  Protocol  (RBP)  III  Score,  relative 
to  the  score  of  a comparable  reference  site. 
Reference  sites  represent  benchmarks  of  the  best 
attainable  conditions  that  can  be  expected  for 
similar  streams  within  a particular  ecological 
region  (Barbour  and  others.  1996).  Thus,  the 
delineation  of  ecological  regions  and  the 
identification  of  similar  streams  and  appropriate 
reference  sites  are  critical  components  of 
bioassessment  methodolosv 

Barbour  and  others  (1996)  describe  the 
process  of  charactenzmg  regional  reference 
conditions  as  consisting  of  three  steps: 


1 


1.  Classification  of  streams  based  on 
geographic  region  and  stream  t\pe; 

2.  Selection  of  the  best  available  sites  within 
a given  class  as  candidate  reference  sites; 
and 

3.  Characterization  of  reference  conditions 
based  on  biological  survey  data. 

These  authors  state  that  the  objectives  of  the  first 
step  m this  process,  the  classification  of  streams 
based  on  geographic  region  and  stream  t\pe,  are 
to  reduce  the  complexity  of  biological  information 
and  improve  the  resolution  or  sensitivity  of 
biological  surveys  by  partitioning  or  accounting 
for  variation  among  sites.  Furthermore,  they 
suggest  that  “the  best  classification  variables  are 
those  that  are  readily  obtained  from  maps  or 
regional  water  characteristics  such  as  ecoregion, 
gradient,  alkalinity,  and  hardness.”  Stream 
charactenstics  that  are  readily  affected  by  human 
activities  or  occur  as  a biological  response  to 
physical  or  chemical  conditions  (i.e.,  land  use, 
habitat  condition,  or  nutrient  concentrations) 
should  not  be  used  as  classification  variables. 

The  methodology  developed  by  Barbour  and 
others  (1996)  indicates  that  stream  classification 
should  include  variables  that  describe  channel 
features  at  a variety  of  different  spatial  scales. 
Gregory  and  others  (1991)  descnbe  a spatial 
hierarchy  of  dramage  basins  that  consists  of 
channel  features  from  five  spatial  scales: 
(1)  network;  (2)  segment;  (3)  reach;  (4)  channel 
unit,  and  (5)  charmel  subunit  (Table  1).  However, 
the  actual  number  and  definition  of  spatial  scales 
utilized  in  stream  classification  programs  vary  and 
may  be  dictated  by  the  specific  objectives  of  the 
program  and  budgetary  and  manpower 
constraints.  For  example,  the  U.S.  Geological 
Survey's  National  Water-Quality  Assessment 
(NAWQA)  Program  utilizes  a hierarchical  stream 
habitat  assessment  design  that  includes 
charactenzing  habitat  at  four  spatial  scales; 
(1)  basin,  (2)  stream  segment;  (3)  stream  reach; 
and  (4)  microhabitat. 

Currently,  the  Pa.  DEP  classifies  streams 
based  on  ecoregion  (or  subregion)  and  stream 
order  when  conducting  bioassessments  and 
ecological  significance  surveys  of  special 
protection  candidate  streams  (Shertzer,  1992). 


The  classification  and  sampling  techniques  used 
by  the  Pa.  DEP  effectively  incorporate 
classification  variables  at  the  segment  and  channel 
unit  spatial  scales,  as  described  in  Table  1,  by 
including  ecoregion  designation  and  by  limiting 
biological  data  collection  to  riffles.  However, 
stream  order,  which  is  used  to  characterize 
streams  at  the  reach  spatial  scale,  does  not 
characterize  stream  reach  morphology,  which  is  a 
critical  component  of  stream  ecosystems.  Stream 
reach  morphology  influences:  (1)  instream  habitat 
conditions  and  hydrograph  characteristics;  (2)  the 
fluvial  development  of  valley  floor  landforms 
such  as  flood  plains,  natural  levees,  and  terraces; 
(3)  riparian  vegetative  communities  and  the 
quantity  and  quality  of  organic  matter  entering 
streams;  and  ultimately,  (4)  aquatic  biological 
community  structure.  For  reviews  of  the 
influence  of  morphological  conditions  on  stream 
and  riparian  ecosystems,  see  Varmote  and  others 
(1980)  and  Gregory'  and  others  (1991). 

The  Rosgen  Stream  Classification  System 
(Rosgen,  1996),  which  is  commonly  used  by  U.S. 
Department  of  Agriculture  (USDA)  Forest 
Service  hydrologists,  provides  a means  by  which 
morphologically-similar  stream  reaches  can  be 
identified.  This  system  utilizes  stream  channel 
features  that  allow  for  the  classification  and 
assessment  of  stream  reach  morphology  at  four 
different  levels  of  resolution.  Rosgen  Level  I and 
II  inventories  are  designed  for  classifying  stream 
reaches.  Level  I classification  is  a coarse-scale 
geomorphic  characterization  that  can  be 
conducted  using  information  taken  from 
topographic  and  landform  maps  and  from  aerial 
photographs.  Level  II  classification,  based  on 
field  data,  is  a more  detailed  morphological 
description  that  provides  a quantitative 
classification  of  stream  reaches.  Rosgen  Level  I 
and  II  inventories  are  summarized  in  Table  2. 

Rosgen  Level  III  and  IV  inventories  are  very 
detailed  and  build  upon  the  information  gathered 
during  Level  Til  inventones.  Level  III 
inventories  are  designed  for  assessing  the  existing 
condition  or  the  “state”  of  specific  stream  reaches. 
Level  IV  inventories  are  used  for  establishing 
empirical  relationships  pertaining  to  the  prediction 
of  Manning’s  “n”  values  and  studying  sediment 
transport  processes  and  hydraulic  geometry 
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Table  1.  Spatial  Scales  Included  in  the  Hierarchical  Organization  of  Drainage  Basins  (Described  by 
Gregory  and  others.  1991) 


Spatial  Scale 

Approximate 
Number  of 
Channel  Widths  in 
Length 

Delineation  Features 

Network 

100,000 

Network  of  fluvial  and  valley  floor  landforms  within  a 
basin  e.xtending  from  headwaters  to  large  rivers  or 
estuaries 

Segment 

1,000  - 10.000 

Regional  landform  patterns  diat  reflect  different  agents  of 
landscape  formation  (scale  of  10  to  100  Km),  major 
topographic  discontinuities  such  as  high  gradient 
mountain  rivers  or  low  gradient  lowland  rivers 

Reach 

10  - 100 

T\pe  and  degree  of  local  constraint  imposed  on  the 
channel  and  valley  floor  by  geomorphic  features 

Channel  Unit 

1 

Water-surface  slope,  channel  width/depth  ratio,  extent  of 
turbulent,  high-velocity  flow;  Upically  classified  as 
riffles,  pools,  and  runs 

Channel  Subunits 

0.1 

Local  hydrauhc  featmes,  at  scales  less  than  the  active 
channel  width,  that  are  created  by  boulders,  logs,  or  gravel 
bars  and  exist  at  flows  less  than  those  required  to  fill  the 
active  channel 

Table  2.  Summary  of  Rosgen  Level  I and  II  Inventories 


Rosgen  Level  I (Geomorphic  Characterization) 

Primary  Information  Sources: 

Topographic  and  landform  maps 

Aerial  photography 

Delineation  Criteria: 

Channel  slope  (valley  slope/sinuosityO 

Charmel  shape  (narrow/deep  vs.  wide/shallow) 

Channel  patterns  (single  thread,  multiple  thread,  anastamosed) 

Stream  Reach  Classes: 

Stream  t>pes  "A"  through  "G" 

Rosgen  Level  II  (Morphological  Description) 

Primary  Information  Sources: 

Field  Data 

Delineation  Criteria: 

Entrenchment  ratio 

Width/depth  ratio 

Sinuosit\' 

Channel  slope 

Charmel  matenals 

Stream  Reach  Classes: 

Stream  t>pes  "Ala...A6a"  through  "Glc...G6c" 

of  streams  based  on  stream  type  (as  determined 
from  Level  II  data)  and  stream  state  (based  on 
Level  III  data).  These  relationships  can  then  be 
extrapolated  to  similar  stream  reaches  for  which 
Level  IV  data  are  not  available. 

To  accurately  classify  candidate  and  reference 
streams  for  the  purpose  of  comparing  their 
macroinvertebrate  communities,  it  is  essential  that 
these  streams  be  similar  at  a variety  of  spatial 
scales.  The  author  hypothesizes  that  the  Pa.  DEP 
could  strengthen  the  procedure  that  is  currently 
used  for  identifying  comparable  Special 
Protection  candidate  and  reference  sites  by 
incorporating  Rosgen  Level  I or  II  stream  reach 
classification  designation  into  this  procedure.  The 
primary  objective  of  this  project  is  to  evaluate  the 
Rosgen  classification  system  as  a method  for 
partitioning  variation  among  the  benthic 
macroinvertebrate  communities  and  for 
characterizing  HQ-CWFs  at  the  stream  reach 
spatial  scale. 

METHODS 

Field  and  Laboratory  Methods 

A total  of  35  sample  sites  located  on 
HQ-CWFs  were  sampled  between  September  1 1 
and  November  11,  1996  (Figure  1 and  Table  3). 
Sample  sites  were  located  in  the  Pennsylvania 
portion  of  the  Susquehanna  River  Basin  and 
represented  6 ecoregions  and  15  subregions 
(Figure  2).  At  each  sample  site,  stream  reach 
morphology  and  chemical  water  quality 
conditions  were  documented  in  the  field,  and 
benthic  macroinvertebrate  samples  were  collected 
for  analysis  in  the  laboratory. 

Chemical  water  quality 

Field  water  quality  measurements  included 
water  temperature,  dissolved  oxygen, 
conductivity,  pH,  alkalinity,  and  acidity. 
Dissolved  oxygen  was  measured  using  a YSI 
Model  55  dissolved  oxygen  meter.  Conductivity 
and  pH  were  measured  using  Cole-Parmer 
Model  1481-61  and  Model  5996-70  meters. 


respectively.  Alkalinity  was  measured  by  titrating 
a known  volume  of  sample  water  to  pH  4.5  with 
0.02N  H2SO4.  Acidity  was  measured  by  titrating 
a known  volume  of  sample  water  to  pH  8.3  with 
0.02N  NaOH. 

Benthic  macroinvertebrate  sampling 

Macroinvertebrate  communities  were  sampled 
using  a 0.3-m  D-ffame  kick  net  to  collect 
organisms  dislodged  from  riffle  areas  by  physical 
agitation  of  approximately  0.3  m^  of  the 
streambed.  Two  samples  were  collected  at  each 
site,  one  from  an  area  of  high  velocity  and  the 
other  from  an  area  of  lower  velocity.  The  two 
samples  were  composited  and  preserved  in  a 
solution  of  isopropyl  alcohol  and  glycerin  for 
laboratory  analysis.  In  the  laboratory, 
100-organism  subsamples  were  taken  from 
composite  samples  using  a gridded  pan  and  a 
random  numbers  table.  The  organisms  contained 
in  the  subsamples  were  identified  to  genus  (except 
for  Chironomidae  and  Simuliidae)  using 
taxonomic  keys  by  Merrit  and  Cummins  (1996) 
and  Peckarsky  and  others  (1990)  and  enumerated. 

Stream  reach  morphology 

At  each  sample  site  (sample  reach),  stream 
reach  morphological  conditions  were  character- 
ized using  Rosgen  Level  II  classification  data 
(Rosgen,  1996).  Level  II  classification  data 
consisted  of  field  measurements  of  channel  cross- 
sectional  dimensions,  sinuosity,  slope,  and  the 
size  of  channel  materials.  Channel  cross-sectional 
dimensions,  sinuosity,  and  slope  were  determined 
using  a surveying  rod,  level,  and  tape  measure.  At 
each  site,  sinuosity,  slope,  and  channel  materials 
were  assessed  over  the  entire  sample  reach. 
Sample  reaches  were  approximately  equal  to  20 
bankflill  channel  widths  in  length.  Sample  reach 
end-points  were  delineated  such  that  reaches 
started  and  ended  in  the  same  type  of  charmel  unit 
(e.g.,  head-of-riffle  to  head-of-riffle). 

At  each  sample  site,  a stream  channel  cross 
section,  located  in  a relatively  straight, 
representative  segment  of  the  sample  reach,  was 
surveyed.  While  surveying  cross  sections,  field 
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Figure  1.  Locations  of  Sanple  Sites  in  the  Susquehanna  River  Basin 


Table  3.  Descriptions  of  35  Sample  Site  Locations 
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Figure  2.  Ecoregiom  of  Pennsylvania  and  the  Locations  of  Sanple  Sites 


indicators  of  bankfull  stage  such  as  those 
described  by  Dunne  and  Leopold  (1978), 
Harrelson  and  others  (1994),  and  Rosgen  (1996) 
were  recorded.  Some  of  these  indicators  included; 
floodplain  elevations,  the  tops  of  point  bars, 
breaks  in  streambank  slope,  changes  in 
streambank  particle  size,  and  staining  of  rocks. 
Bankfull  channel  dimensions  and  estimated 
bankfull  discharge  were  determined  from 
surveyed  cross-sectional  and  surface  water  slope 
data  and  channel  substrate  composition 
information  using  XSPRO  stream  channel  cross- 
sectional  analysis  software  (Grant  and  others, 
1992).  The  Thome  and  Zevenbergen  model  in 
XSPRO  was  used  to  estimate  stream  discharge. 
The  model  was  mn  using  the  d84  particle  size  from 
pebble-count  data  and  sample  reach  surface  water 
slope  (described  below).  At  sample  sites  that 
were  substantially  influenced  by  bedrock,  the  d84 
values  used  in  the  Thorne  and  Zevenbergen  model 
were  adjusted  until  reasonable  Manning's  ‘ft” 
values  (based  on  Figure  8-2  of  Rosgen  (1996)) 
were  obtained.  The  bankftill  stage  elevation  of 
each  sample  site  was  determined  in  the  office, 
based  on  a review  of  the  following:  (1)  channel 
cross-sectional  survey  data  and  field  notes;  (2)  a 
comparison  of  computer-generated  printouts  of 
estimated  discharge  at  0.05-foot  stage  increments 
to  average  values  of  bankfull  channel  discharge  as 
a function  of  drainage  area  for  streams  in 
Pennsylvania  reported  in  Dunne  and  Leopold 
(1978);  and  (3)  photographs  taken  at  sample  sites. 

Entrenchment  ratio  describes  the  vertical 
containment  of  the  stream  channel  (the  degree  to 
which  it  is  incised  in  the  valley  floor)  and 
provides  insight  into  how  stream  energy  is 
dissipated  during  high  flows.  For  example,  the 
top-of-bank  elevations  of  entrenched  streams 
(stream  types  F and  G)  are  much  higher  than 
bankfull  stage  (Figure  3).  At  flows  greater  than 
bankfull  stage,  streamflow  is  contained  within  the 
streambanks,  and  stream  depth  increases  much 
faster  than  width.  During  floods,  this  containment 
of  streamflow,  which  results  in  stream  energy 
being  exerted  on  both  streambanks,  is  partly 
responsible  for  the  high  bank  erosion  rates 
observed  in  entrenched  streams  (Rosgen.  1996). 
In  contrast  to  entrenched  streams,  slightly 
entrenched  streams  (stream  types  C and  E) 


overtop  their  streambanks  at  flows  greater  than 
bankfull  stage,  and  stream  width  increases 
substantially  and  stream  energy  is  dissipated  as 
streamflow  moves  onto  the  stream’s  relatively 
broad  flood  plain.  The  vertical  containment  of 
moderately  entrenched  streams  (type  B streams)  is 
intermediate  to  that  of  entrenched  and  slightly 
entrenched  streams. 

Entrenchment  and  width/depth  ratios  were 
derived  from  channel  cross-sectional  data. 
Entrenchment  ratio  was  calculated  as  the  ratio  of 
the  width  of  the  flood-prone  area  to  the  width  of 
the  bankfull  channel.  The  flood-prone  area  width 
was  measured  at  the  elevation  that  corresponded 
to  twice  the  maximum  depth  of  the  bankfull 
channel  (Figure  4).  Width/depth  ratio  was 
calculated  as  the  ratio  of  the  surface  width  to  the 
mean  depth  of  the  bankfull  channel. 

Sinuosity  was  calculated  as  the  ratio  of  stream 
channel  length  to  valley  length  of  the  study  reach. 
Stream  channel  length  was  determined  in  the  field 
by  measuring  the  length  of  each  channel  unit 
(riffle,  run,  or  pool)  in  the  study  reach  using  a tape 
measure  extended  in  alignment  with  the  channel 
thalweg  (the  deepest  part  of  the  channel).  Valley 
length  was  determined  by  measuring  the  straight- 
line  distance  between  the  upstream  and 
downstream  end-points  of  the  reach,  which  were 
surveyed  in  the  field,  and  plotted  in  the  office. 
Water  surface  slope  was  calculated  as  the  change 
in  surface  water  elevation  per  unit  stream  length. 
Surface  water  elevation  at  the  upstream  and 
downstream  end-points  of  each  study  reach  were 
surveyed  in  the  field.  Stream  length  was 
determined  in  the  field  using  a tape  measure,  as 
described  above. 

Dominant  channel  materials  were  determined 
based  on  the  particle  size  distribution  of  bankfull 
channel  materials  using  the  modified  Wolman 
“pebble-count”method,  as  described  by  Rosgen 
(1996).  The  intermediate  axis  diameter  of 
approximately  100  randomly-selected  substrate 
particles  was  measured  and  tabulated  in  the  field. 
The  number  of  particles  selected  from  a given 
channel  unit  (riffle,  run,  or  pool)  was  based  on  the 
channel  unit  configuration  of  the  stream  reach, 
such  that  channel  units  were  sampled  in 
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I'i^iirc  jf.  Examples  of  Entrenched,  Moderately  Entrenched,  and  Slii’htly  Entrenched  Stream  Types 


Fisiire  4.  Illustration  of  Entrenchment  Ratio  Calculated  From  Bankfull  and  Flood-Prone  Channel 
Widths 
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proportion  to  their  occurrence  in  the  reach.  The 
dominant  particle  size  of  the  stream  reach  was 
determined  as  the  indexed  dso  value  from  a 
cumulative  percent  plot  of  the  pebble-count  data 
(Appendix  A).  Rosgen  Level  II  stream  t\pe 
designations  were  determined  using  the  cntena 
summanzed  m Table  4. 

Ecoreqion  designation,  stream  order,  and 

watershed  characteristics 


Sample  site  ecoregion  (subregion)  designation 
and  watershed  charactenstics  (e.g.,  dramage  area, 
percent  ecoregion  composition)  were  determined 
using  PC  Arc  View  spatial  analyst  geographic 
information  system  (GIS)  software  and  readily 
available  geographic  data.  Stream  order  (Strahler 
method)  was  delineated  from  U.S.  Geological 
Survey  7.5  minute  quadrangle  maps.  Intermittent 
streams  were  not  included  in  stream  order 
designations. 

Data  Analysis 

All  statistical  tests  and  simple  linear 
correlation  analyses  were  conducted  in 
accordance  with  Zar  (1984)  using  MINITAB 
statistical  software  (MINITAB,  1996).  Box-and- 
vvhisker  plots  were  generated  using  this  same 
software.  Box-and-whisker  outliers  are  points 
outside  of  the  lower  and  upper  limits  of  the 
sample  group.  Significant  results  were  reported  at 
the  95  percent  confidence  level  (p<0.05).  When 
statistically-significant  analysis  of  variance 
(ANOVA)  results  w^ere  obtained,  differences 
among  the  sample  groups  were  evaluated  using 
Fisher's  pairwise  multiple  comparisons  and  box- 
and-whisker  plots.  Box-and-whisker  plots 
graphically  display  the  median,  minimum, 
maximum.  25^  and  75‘^  percentiles  (also  referred 
to  as  the  V''  and  3'^^^  quartiles,  respectively),  and 
outliers  of  a given  sample  group.  Sample  group 
limits  were  determined  as  follows; 

Lower  Limit:  quartile  - 1.5 

(3'^'^  quartile  - quartile) 

Upper  Limit:  3'^'^  quartile  - 1.5 

(3'^‘^  quartile  - U*  quartile) 


Some  of  the  statistical  procedures  used  in  the 
study  required  grouping  sample  sites  into  drainage 
area  and  surface  water  slope  classes.  These 
groupings  were  based  on  the  sites'  dramage  area 
and  surface  water  slope  values  relative  to  the  other 
sample  sites  in  the  data  set.  For  example,  sample 
sites  were  grouped  into  one  of  four  dramage  area 
classes.  These  classes  were  based  on  the 
following  critena: 


Dramage 

Area  Class  Cntena 


1 Dramage  area  size  < the 

U'  quartile  value  of  the  data 
set 

2 Dramage  area  size  > or 

equal  to  the  U'  quartile 

value  and  < the  median  of 
the  data  set 

3 Dramage  area  size  > or 

equal  to  the  median  and  < 
the  3’’'^  quartile  value  of  the 
data  set 

4 Dramage  area  size  > or 

equal  to  the  3'''*  quartile 

value  of  the  data  set 

A similar  procedure  was  used  to  group  sample 
sites  into  surface  water  slope  classes. 


Detrended  correspondence  analysis  (DCA) 
was  used  to  reveal  outliers  in  the  data  set  based  on 
macromvertebrate  taxonomic  abundance.  Gauch 
(1982)  describes  outliers  as  peculiar  samples 
whose  relationships  to  other  samples  in  the  data 
set  are  not  expressed  by  the  information  in  the 
data.  Furthermore.  Gauch  states  that,  in  general, 
outliers  tend  to  compress  the  other  samples  m the 
data  set  into  an  undesirably  small  area  of  the 
ordination  and  obscure  the  analysis  of  the  data  set 
as  a whole,  and  thus,  should  be  deleted  from  the 
data  set.  Since  the  focus  of  this  study  is  on  the 
relationships  between  stream  reach  morpholog\’ 
and  the  benthic  macroinvertebrate  communities  of 
HQ-CWFs.  samples  that  were  pnmanly 
influenced  by  water  quality  conditions  or  that  had 
obviously  degraded  channel  morphologv  or 
macroinvertebrate  communities  were  removed 
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Values  of  Entrenchment  and  Sinuosity  can  vary  by  +/-0.2  units;  values  for  Uldth/Depth  ratios  can  vary  by+/-  2.0  units.  Source:  Modified  troin  Rosgen  (1996). 


from  the  analysis.  DCA  ordinations  were 
produced  using  software  developed  by  Kovach 
(1993).  Rare  taxa.  those  present  at  less  than 
5 percent  of  the  sites,  were  removed  pnor  to 
running  ordination  analyses  (Gauch.  1982). 

Five  metrics  were  used  to  descnbe  sample  site 
benthic  macromvertebrate  communities.  These 


metrics,  which  included  taxa  nchness.  modified 
EPT  index,  modified  Hilsenhoff  index,  percent 
dominant  ta.xa.  and  modified  percent  ma\flv 
index,  are  summanzed  in  Table  5.  The  Pa.  DEP 
was  considering  this  suite  of  metncs  for  assessing 
the  ecological  significance  of  Special  Protection 
candidate-streams  at  the  start  of  this  project 
(Pa.  DEP.  Unpublished). 


Table  5.  Summary  of  the  Bioassessment  Metrics  Used  to  Describe  Sample  Site  Benthic 
Macroinvertebrate  Communities 


Metric 

Description 

Taxa  Richness 

Total  number  of  different  taxa. 

Modified  EPT  Index 

Number  of  pollution  intolerant  (Hilsenhoff  tolerance  value  < 5) 
Ephemeroptera  (ma\fly),  Plecoptera  (stonefly),  and  Trichoptera  (caddisfly) 
taxa. 

Modified  Hilsenhoff  Biotic  Index 

Hilsenhoff  Biotic  Index  of  pollution  tolerance  with  pollution  tolerance 
values  modified  for  Pennsjivania. 

Percent  Dominant  Taxa 

Relative  abundance  of  the  dominant  taxa  expressed  as  percent  of  the  100- 
organism  subsample. 

Modified  Percent  Mavdly  Index 

Relative  abundance  of  pollution  intolerant  (Hilsenhoff  tolerance  value  < 5) 
mavdlies  expressed  as  percent  of  the  100-organism  subsample. 

RESULTS  AND  DISCUSSION 
Stream  Classification 

Rosgen  Level  I classification  delineates  the 
35  sample  sites  included  in  this  study  into  four 
different  stream  t\pes.  The  majority  of  the 
sample  sites  are  either  B or  C stream  t>pes,  which 
together  account  for  89  percent  of  the  sites 
included  m the  study.  Rosgen  Level  II 
classification  refines  the  major  stream  t\pes  (t\pes 
A through  G)  into  subcategories  based  on  channel 
materials  and  water  surface  slope.  For  example,  a 
Level  II  B3c  stream  reach  is  a B stream  type  with 
cobble  substrate  and  a surface  water  slope  of  less 
than  0.020  (see  Table  4).  Rosgen  Level  II  stream 
classification  further  delineates  the  35  sample 
sites  into  1 1 different  stream  t\pes  consisting 
pnmarily  of  types  C4  (37  percent),  C3 
(17  percent),  C4b  (11  percent),  and  B3  (9  percent) 
(Figure  5). 

Due  to  the  excessive  fragmentation  of  the  data 
set  and  the  low  number  of  streams  in  many  of  the 
Level  II  stream  t\pe  categories,  a tliird  level  of 
classification  was  developed.  This  additional 
level  of  classification,  referred  to  as  the  “modified 
Rosgen”  level,  produced  stream  type  designations 
intermediate  to  those  of  Rosgen  Levels  I and  II. 
Modified  Rosgen  classification  involved 
classifying  sample  sites  based  on  their  Rosgen 
Level  II  designation  without  the  surface  water 
slope  subcategory.  Modified  Rosgen  classifi- 
cation partitioned  the  35  sample  sites  into  one  of 
eight  stream  fype  categories  (Figure  5).  Rosgen 
stream  reach  classification  information  is 
summanzed  m Tables  6 and  7 and  Appendix  A. 
Relationships  between  estimated  bankfull 
discharge  and  drainage  area  and  between  bankfull 
channel  dimensions  and  drainage  area  of  sample 
sites  are  shown  in  Figures  6 and  7,  respectively. 
Water  quality  data  are  listed  in  Appendix  B. 

Benthic  Macroinvertebrate  Community 
Structure 

A total  of  83  taxa  from  43  families  were 
included  m the  macroinvertebrate  data  collected  at 
the  35  sample  sites  (Appendix  C).  Figure  8 shows 
the  DC  A ordination  of  35  sample  sites  (data  set  1) 
based  on  macroinvertebrate  taxa  abundance.  This 


ordination  suggests  LTRT,  HAMU,  HAML, 
WILL,  HORS,  HUl'I,  and  HARV  (see  Table  3) 
may  be  outliers  in  the  data  set.  Correlation 
analysis  of  channel  morphology  and  water  quality 
data  with  sample  site  DCA  ordination  scores 
suggest  that  macroinvertebrate  community 
structure  of  data  set  1 is  primarily  influenced  by 
conductivity  (r  = 0.775,  p < 0.001)  and  alkalinity 
(r  = 0.737,  p < 0.001),  and  to  a lesser  extent, 
degraded  channel  substrate  (Figure  9).  LTRT. 
HAMU,  HAML,  WILL,  HORS,  and  FIUTT  were 
determined  to  be  outliers  based  on  the  rationale 
described  in  Table  8.  HARV  was  retained  since 
it;  (1)  IS  not  influenced  by  limestone  geology; 
(2)  is  characterized  by  stable,  natural  stream  reach 
morphology;  and  (3)  supports  a healthy 
macroinvertebrate  community. 

A second  data  set,  consisting  of  29  sites  (data 
set  2),  was  constructed  by  removing  the  six 
outliers  mentioned  above  from  the  original  of 
35  sample  sites  (data  set  1).  DCA  performed  on 
data  set  2 produced  a DCA  ordination  of 
29  sample  sites  that  does  not  suggest  the  presence 
of  any  additional  outliers  (Figure  10).  Correlation 
analysis  of  DCA  axis  scores,  morphological 
characteristics,  and  water  quality  data  indicate  that 
macroinvertebrate  community  structure  in  data 
set  2 is  predominantly  influenced  by  dissolved 
oxygen,  percent  cobble  substrate,  acidity,  and 
percent  riffle.  The  metric  scores  of  all  35  sample 
sites  are  listed  in  Table  9.  The  metric  scores  of 
data  set  1 (35  sites)  and  data  set  2 (29  sites, 
outliers  identified  by  DCA  removed)  are 
summanzed  in  Table  10  and  Figure  11. 

Relationships  Between  Stream  Type  and 
Metric  Scores 

Twenty-seven  of  the  29  sample  sites  included 
in  data  set  2 are  classified  as  Rosgen  Level  I type 
B or  C streams  (Figure  12).  These  27  sites  (data 
set  3)  were  used  to  evaluate  the  effectiveness  of 
ecoregion,  stream  order,  drainage  area,  and 
Rosgen  stream  type  designation,  predominant 
channel  material,  and  surface  water  slope  as 
stream  classification  variables.  Data  set  3 stream 
classification  and  benthic  macroinvertebrate 
community  data  are  summarized  in  Table  1 1 . The 
results  of  statistical  analyses  to  test  for  significant 
differences  in  metric  scores  based  on  stream 
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NUMBER  OF  SITES  NUMBER  OF  SITES 


B C E F 

ROSGEN  LEVEL  I STREAM  TYPE 


B3a  B3  B3c  B4a  B5c  C3  C4b  C4  E4  E5b  F4 

ROSGEN  LEVEL  II  STREAM  TYPE 


20  -r 


MODIFIED  ROSGEN  STREAM  TYPE 


Figure  5.  Histograms  of  35  Sample  Sites  Based  on  Rosgen  Level  I,  Level  //,  and  Modified  Rosgen 
Stream  Type  Designation 
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Table  6.  Summary  of  Rosgen  Stream  Reach  Classification  Data  of  35  Sample  Sites 


Sample 

Site 

Entrenchment 

WD  Ratio 

Sinuosity 

Slope 

Channel 

Material 

(dso) 

Rosgen 
Level  1 

Class 

Rosgen 
Level  II 

Class 

Modified 

Rosgen 

Class 

ASPH 

1.8 

27.5 

1.1 

0.020 

Cobble 

B 

B3 

B3 

BALD 

1.5 

22.2 

1.1 

0.029 

Cobble 

B 

B3 

B3 

CORK 

>2.2 

24.0 

1.2 

0.019 

Gravel 

C 

C4 

C4 

DEEP 

2.0 

22.6 

1.1 

0.036 

Cobble 

B 

B3 

B3 

DWMN 

>2.2 

24.3 

1.4 

0.029 

Gravel 

C 

C4b 

C4 

EBCC 

>2.2 

14.1 

1.0 

0.002 

Gravel 

C 

C4 

C4 

EBLR 

>2.2 

49.4 

1,0 

0.010 

Cobble 

c 

C3 

C3 

EBRB 

1.5 

17.0 

1.0 

0.019 

Cobble 

B 

B3c 

B3 

GAZZ 

>2.2 

14.2 

1.4 

0.008 

Gravel 

C 

C4 

C4 

GRLE 

>2.2 

14.2 

1,4 

0.014 

Cobble 

c 

C3 

C3 

GRNL 

>2.2 

14.8 

1.5 

0.014 

Gravel 

c 

C4 

C4 

GRNU 

>2.2 

14.6 

1.1 

0.012 

Gravel 

c 

C4 

C4 

HAML 

1.3 

19.9 

1.1 

0.003 

Gravel 

F 

F4 

F4 

HAiVrU 

1.9 

10.3 

1.0 

0.003 

Sand 

B 

B5c 

B5 

HARV 

>2.2 

30.1 

1.2 

0.005 

Gravel 

C 

C4 

C4 

HORS 

>2.2 

7.7 

1.2 

0.025 

Sand 

E 

E5b 

E5 

HUTT 

1.8 

28.3 

1.1 

0.040 

Gravel 

B 

B4a 

B4 

IRON 

>2.2 

24.6 

1.1 

0,019 

Gravel 

C 

C4 

C4 

LEBO 

>2,2 

22.5 

1.2 

0.015 

Cobble 

C 

C3 

C3 

LMOR 

>2.2 

21.2 

1.2 

0.027 

Gravel 

C 

C4b 

C4 

LTRT 

1.7 

13.8 

1.1 

0.006 

Sand 

B 

B5c 

B5 

MNTN 

>2.2 

19.2 

1.1 

0.005 

Gravel 

C 

C4 

C4 

MURY 

>2.2 

16,0 

1.0 

0.014 

Gravel 

C 

C4 

C4 

PEQU 

>2.2 

15.3 

1.0 

0.017 

Gravel 

C 

C4 

C4 

RICE 

>2.2 

17.6 

1.2 

0.040 

Gravel 

C 

C4b 

C4 

SIXM 

>2.2 

34.7 

1.0 

0.016 

Cobble 

C 

C3 

C3 

STSN 

>2.2 

12.5 

1.2 

0.009 

Cobble 

C 

C3 

C3 

TIPT 

>2.2 

24.0 

1,1 

0.010 

Cobble 

C 

C3 

C3 

TOML 

>2.2 

8.3 

1.0 

0.001 

Gravel 

E 

E4 

E4 

TOMU 

>2.2 

33.0 

1,1 

0,020 

Gravel 

C 

C4b 

C4 

UEBC 

2.0 

8,0 

1.2 

0.014 

Gravel 

E 

E4 

E4 

UHAR 

2.2 

13.5 

1,1 

0,041 

Cobble 

B 

B3a 

B3 

URRB 

>2.2 

23.7 

1.1 

0,018 

Gravel 

C 

C4 

C4 

WILL 

>2.2 

36.9 

1.2 

0.005 

Gravel 

C 

C4 

C4 

WTDR 

>2.2 

26.8 

1,3 

0.004 

Gravel 

C 

C4 

C4 

16 


Table  7.  Distribution  of  35  Sample  Sites  by  Rosfjen  Level  II  Stream  Type  and  Ecore^ion  Designation 
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64c  Piedmont  Uplands  69a  Forested  hills  and  mountains 

69b  Uplands  and  valleys  of  mixed  land  Inse 


0 1 10  100 


DRAINAGE  AREA  (mi') 


I 8 Sample  Sites Power  (Dunne  and  Leopold  (1978))  I 


Figure  6.  Scatter  Plot  of  Estimated  Bankfull  Discharge  as  a Function  of  Drainage  Area  of  Sample 
Sites  (Solid  line  represents  average  values  for  streams  in  Pennsylvania  based  on  the  above 
equation  reported  by  Dunne  and  Leopold  (1978)) 
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1 10 
DRAINAGE  AREA  (mi') 

(ai 


DRAINAGE  AREA  |mi') 

(bj 


9 Sample  Sites 


- Po^er  (Dunne  and  Leopold  (1979^) 


(c) 


Figure  7.  Scatter  Plot  of  (a)  Bankfull  Channel  Area,  (b)  Bankfull  Channel  Width,  and  (c)  Bankfiill 
Channel  Average  Depth  as  Functions  of  Drainage  Area  of  Sample  Sites  (solid  line 
represents  average  values  for  streams  in  the  Eastern  United  States  estimated  from  graphs  in 
Dunne  and  Leopold  (1978) 
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Figure  8.  Detrended  Correspondence  Analysis  (DCA)  Ordination  of  35  Sample  Sites  Based  on 
Macroinvertebrate  Taxonomic  Abundance 
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Table  8.  Summary  of  the  Six  Outlier  Sample  Sites  and  the  Justification  for  Their  Removal  From  the 
Original  Data  Set  of  35  Sample  Sites 


Site 

Justification  for  Removal  From  Data  Set 

LTRT 

Elevated  conductivity  and  alkalinity,  geochemically  influenced. 

Site  is  located  ia  and  entire  watershed  consists  of.  Ecoregion  67a  (limestone  or  dolomite  valleys 
subregion). 

Degraded  channel  morphology  (G5c  stream  type  in  disequilibrium). 

HAiVIU 

Elevated  conductivity-  and  alkalinity',  geochemically  influenced. 

Site  is  located  in.  and  5 1 percent  of  its  watershed  consists  of,  Ecoregion  67a  (limestone  or  dolomite 
valleys  subregion). 

Degraded  channel  morphologv'  (G5c  stream  type  in  disequilibrium). 

HAML 

Elevated  conductivity'  and  alkalinity,  geochemically  influenced. 

60  percent  of  the  site's  watershed  consists  of  Ecoregion  67a  (Umestone  or  dolomite  valleys 
subregion). 

Degraded  channel  morphology  (F4  stream  type  in  disequilibrium). 

WILL 

Elevated  conductivity  and  alkalinity',  geochemically  influenced. 

12  percent  of  the  site's  watershed  consists  of  Ecoregion  67a  (hmestone  or  dolomite  valleys 
subregion). 

HORS 

Degraded  channel  substrate  composition  (52  percent  sand). 

Poor  metric  scores. 

HUTT 

Poor  metric  scores. 
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Figure  10.  DCA  Ordination  of  Data  Set  2 Sample  Sites  (Data  set  2 consists  of  29  sites  and  was 

constructed  by  removing  the  six  outliers,  identified  by  DCA,  from  the  original  data  set  of 
25  sample  sites  (data  set  1).) 


Table  9.  Metric  Scores  of  35  Sample  Sites 


Site 

Tax  Rich 

Mod  EPT  1 

Mod  Mils 

% Dom 

Mod  % May 

ASPH 

25 

15 

2.37 

21 

21 

BALD 

19 

10  1 

2.28 

16 

15 

CORK 

16 

11 

1.66 

27 

36 

DEEP 

24 

15 

1.52 

25 

25 

DWMN 

7 

4 

2.81 

29 

0 

EBCC 

17 

10 

3.43 

29 

25 

EBLR 

20 

12 

2.48 

19 

44 

EBRB 

28 

16 

1.98 

13 

20 

GAZZ 

25 

9 

3.75 

14 

10 

GRLE 

13 

5 

3.10 

24 

8 

GRNL 

22 

12 

3.00 

18 

15 

GRNU 

20 

10 

2.80 

16 

22 

HAML* 

10 

4 

4.28 

41 

10 

HAMU* 

10 

2 

5.52 

31 

1 

HARV 

16 

10 

3.32 

20 

42 

HORS* 

J 

0 

7.07 

79 

0 

HUTT* 

5 

1.68 

48 

48 

IRON 

21 

10 

2.69 

15 

3 

LEBO 

19 

13 

2.91 

27 

12 

LMOR 

23 

12 

3.56 

16 

18 

LTRT* 

7 

0 

6.54 

56 

0 

MNTN 

21 

12 

3.58 

20 

16 

MURY 

20 

9 

3.00 

21 

17 

PEQU 

19 

8 

3.72 

49 

7 

RICE 

14 

11 

2.07 

22 

15 

SDxM 

17 

9 

2.70 

21 

STSN 

17 

9 

2.38 

29 

2 

TIPT 

21 

11 

2.31 

37 

10 

TOML 

16 

7 

3.41 

12 

0 

TOMU 

17 

10 

2.77 

13 

2 

UEBC 

17 

7 

2.86 

22 

9 

UHAR 

15 

10 

3.08 

21 

23 

URRB 

20 

11 

2.41 

20 

14 

WILL* 

20 

10 

3.23 

27 

45 

WTDR 

22 

13 

2.44 

15 

5 

* Outlier  based  on  ordination  of  DCA  axis  1 and  2 scores 
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Table  10.  Summary  of  Metric  Scores  of  Data  Sets  1 (35  Sample  Sites)  and  2 (29  Sample  Sites, 
6 Outliers  Identified  by  DCA  Removed) 


Taxonomic 

Richness 

Modified 

EPT 

Index 

Modified 

Hilsenhoff 

Index 

Percent 

Dominant 

Index 

Modified 
Percent  Mayfly 
Index 

Data  Set 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

Number  of  Sites 

35 

29 

35 

29 

35 

29 

35 

29 

35 

29 

Minimum 

D 

7 

0 

4 

1.52 

1.52 

12 

12 

0 

0 

Maxunum 

28 

28 

16 

16 

7.07 

3.75 

79 

49 

48 

44 

Mean 

17.3 

19.0 

9.1 

10.4 

3.10 

2.80 

26.1 

21.8 

15.6 

15.2 

Median 

19.0 

19.0 

10.0 

10.0 

2.90 

2.80 

21.0 

21.0 

14.0 

15.0 

PERCENT  DOMINANT 


Taxonomic  Richness  vs.  Data  Set 


Modified  EPT  Index  vs.  Data  Set 
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Figure  1 1.  Summary  of  Metric  Scores  of  Data  Sets  1 (35  Sample  Sites)  and  2 (29  Sample  Sites, 
6 Outliers  Identified  by  DCA  Removed) 
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Figure  12.  Histograms  of  Data  Set  3 Sample  Sites  (27  Sample  Sites)  Based  on  Rosgen  Level  /, 
Level  II,  and  Modified  Rosgen  Stream  Type  Designation 
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Table  II.  Summary  Table  of  Data  Set  3 Stream  Classification  and  Jienthic  Macroinvertebrate  Community  Data 


Modified 

Percent 

Mayfly 

Index 

— 

£ 

ir-\ 

04 

44 

04 

— 

OC 

04 

04 

0^4 

T 

ro 

04 

OC 

bO 

WO 

ro 

04 

04 

fO 

04 

or 

Percent 

Dominant 

Taxa 

Index 

r-i 

so 

27 

»/^ 

29 

29 

Cb 

ro 

24 

OC 

bC 

04 

27 

bC 

20 

04 

Cb 

or 

22 

04 

29 

r- 

ro 

ro 

04 

20 

yr\ 

Modified 

Hilsenhoff 

Biotic 

Index 

2.37 

2.28 

1 .66 

1.52 

OC 

ro 

-T 

ro 

cc 

TT 

oi 

1.98 

3.75 

3.10 

3.00 

2.80 

04 

r*“, 

69  t 

2.91 

3.56 

3.58 

3.00 

3.72 

2.07 

2.70 

2.38 

2.31 

LLZ 

3.08 

2.41 

2.44 

Modified 
EPT  Index 

*r\ 

01 

— 

'n 

04 

bO 

Cb 

to 

04 

10 

— 

10 

r^. 

04 

04 

Cb 

OC 

— 

C 

C 

- 

o 

o 

— 

ro 

Taxonomic 

Richness 

25 

C^ 

'sO 

24 

r- 

20 

OC 

04 

25 

ro 

22 

04 

bC 

04 

Cb 

0“, 

04 

04 

(N 

c 

or 

r- 

04 

'/O 

CN 

22 

Predominant 

Channel 

Material 

Cobble 

Cobble 

Ij 

> 

Cobble 

Gravel 

Gravel 

Cobble 

Cobble 

Gravel 

Cobble 

Gravel 

Gravel 

13 

> 

ro 

u- 

o 

Gravel 

Cobble 

Gravel 

Gravel 

Gravel 

Gravel 

Gravel 

Cobble 

Cobble 

Cobble 

Gravel 

Cobble 

Gravel 

Gravel 

Surface 

Water 

Slope 

Class^ 

TT 

— 

04 

— 

04 

04 

04 

— 

04 

•c 

— 

04 

r^i 

or 

— 

04 

or 

or 

m 

— 

Modified 

Rosgen 

Stream 

Type 

rr\ 

cc 

C4 

B3 

C4 

C4 

C3 

B3 

C4 

C3 

C4 

U 

u 

C4 

ro 

U 

U 

C4 

C4 

U 

C4 

ro 

u 

C3 

C3 

C4 

B3 

C4 

C4 

Rosgen 
Level  II 
Stream 
Type 

r^, 

C4 

C4b 

C4 

ro 

U 

O 

ro 

C4 

C3 

C4 

C4 

C4 

C4 

C3 

C4b 

C4 

C4 

C4 

C4b 

C3 

C3 

C3 

C4b 

B3a 

or 

CJ 

C4 

Rosgen 
Level  1 
Stream 
Type 

/-V-V 

u 

c 

u 

CJ 

u 

c 

u 

u 

u 

CJ 

O 

u 

u 

CJ 

a 

o 

u 

CJ 

O 

u 

U 

u 

e 

CJ 

u 

! 

Drainage 

Area 

Class^ 

— 

CO 

-r 

ro 

04 

04 

— 

— 

ro 

04 

or 

ro 

04 

-- 

or 

or 

— 

-- 

04 

CO 

Stream 

Order 

rn 

— 

— 

<N 

r*, 

04 

04 

— 

04 

04 

04 

04 

ro 

04 

ro 

— 

04 

— 

m 

r^. 

04 

- 

- 

04 

04 

Subregion 

62c 

Pc9 

6()b 

62c 

66b 

60b 

62b 

69b 

r- 

bC 

69b 

69b 

60a 

66a 

P69 

62d 

66b 

67d 

64c 

62c 

69a 

67b 

67d 

66b 

60a 

62b 

67c 

Ecoregion 

1 62 

62 

62 

66 

64 

62 

69 

67 

Cb 

bC 

Cb 

bC 

bC 

99 

62 

62 

99 

67 

64 

62 

69 

67 

67 

99 

60 

62 

67 

Sample 

Site 

< 

< 

CORK 

t: 

te 

C 

u 

u 

CC 

UlJ 

GAZZ 

GRNL 

GRNIJ 

IIARV 

IRON 

LEBO 

LMOR 

;2 

§ 

MURY 

O' 

ea 

1^ 

IRICE 

2 

a 

OC 

2: 

OC 

c— . 

Ea 

2 

o 

rs/ 

§ 

5 

cC 

C4 

C 

I— ( 
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Drainage  Area  Classes;  1 less  than  1.99  inP  2 Surface  Water  Slope  Classes;  1 less  than  0.010 

2 1.99  to  3.47  mi^  2 0.010  to  0.015 

3 3,48  to  12.92  mi-  3 0.016  to  0.019 

4 12.93  to  18.54  mi-  4 0.020  to  0,041 


classification  vanables  are  summanzcd  in 
Table  12. 

ANOVA  results  indicate  that  the  percent 
dominant  taxa  and  modified  percent  ma>fl\-  index 
scores  of  data  set  3 sample  sites  differ 
significantly  among  ecoregions  (Figure  13  and 
Table  12).  No  significant  differences  were 
observed  in  any  of  the  five  metrics  included  in  the 
study  based  on  stream  order,  drainage  area  class, 
or  Rosgen  Level  I stream  t\pe  designation 
(Table  12).  ANOVA  results  of  Rosgen  Level  II 
classification  data  show  that  taxa  nchness  scores 
of  B3  and  B3c  streams  are  significantly  higher 
than  those  of  C4b  streams  (Figure  14  (a)).  Rosgen 
Level  II  data  also  indicate  that  modified  percent 
mayfly  index  scores  of  C4b  streams  in  Group  E 
(data  set  3 sites  in  Ecoregions  64  66,  67,  and  69; 
Table  12)  are  significantly  lower  that  those  of  C4 
streams  in  this  group  (Figure  14(b)  and  Table  12). 
However,  it  should  be  noted  that  two  of  the  four 
C4b  sample  sites  (TOMU  and  DWiVIN)  are 
located  m Ecoregion  66b  and  have  the  lowest  field 
alkalinity  (2  mg/L  as  CaCOs)  and  the  highest 
acidity  (14  and  16  mg/L  as  CaCOs,  respectively) 
observed  in  the  study.  Furthermore,  TOMU  and 
DWIVIN  are  the  only  C4b  streams  in  Group  E,  and 
the  two  remaining  C4b  streams  (LMOR  and 
RICE)  are  located  in  Ecoregion  62  and  have 
substantially  higher  modified  percent  mavfly 
index  scores  than  TOMU  and  DWMN.  Thus,  the 
differences  observed  in  the  taxa  richness  and 
modified  percent  mavfly  index  scores  of  Rosgen 
Level  II  stream  tvpes  are  more  readily  attnbuted 
to  ecoregion,  rather  than  stream  reach  channel 
morphology  conditions. 

As  mentioned  previously,  excluding  the 
surface  water  slope  subcategoiy-  designation  and 
classifying  sample  sites  based  on  Level  I stream 
type  and  predominant  channel  material  was  used 
to  classify  streams  into  "modified  Rosgen"  stream 
types.  Data  set  3 ANOVA  results  indicate  that  the 
modified  EPT  mdex  scores  of  modified  Rosgen 
B3  streams  are  significantly  higher  than  C3  and 
C4  streams  (Figure  15(a)  and  Table  12).  In 
addition,  the  percent  dominant  taxa  index  scores 
of  Group  B modified  Rosgen  C3  streams  are 
significantly  higher  than  those  of  B3  and  C4 
streams  in  this  group  (Figure  15(b)  and  Table  12). 


To  further  investigate  the  influence  of  stream 
reach  surface  water  slope  and  channel  substrate 
composition  on  benthic  macroinvertebrate 
community  structure  and  composition,  statistical 
analyses  were  run  on  data  set  3 sample  sites  based 
on  surface  water  slope  class  and  predominant 
channel  matenal  category'.  ANOVA  results 
indicate  that  there  is  no  significant  difference  in 
any  of  the  five  metrics  included  in  the  study  based 
on  surface  water  slope  class  (Table  12).  T-test 
results  indicate  the  mean  modified  Hilsenhoff 
index  score  of  gravel  streams  is  significantly 
higher  than  that  of  cobble  streams  (Table  12  and 
Figure  16).  However,  all  of  the  sites  included  in 
this  analysis  tended  to  have  relatively  low 
modified  Hilsenhoff  index  scores 
(range  = 1 .52  - 3.75),  and  the  ecological 
significance  of  the  observed  difference  in  the 
Hilsenhoff  index  scores  of  cobble  vs.  gravel 
streams  is  negligible. 

Of  the  five  stream  reach  classification 
vanables  (stream  order,  dramage  area.  Rosgen 
Level  I,  Rosgen  Level  11,  and  modified  Rosgen 
stream  type  designation)  included  in  this  study, 
the  modified  Rosgen  stream  type  designation  is 
the  most  effective  classification  variable  for 
partitioning  variation  among  the  benthic 
macromvertebrate  communities  of  HQ-CWTs. 
This  conclusion  is  based  on  the  following; 

1.  No  significant  differences  were  observed 
m any  of  the  five  metncs  included  in  the 
study  based  on  stream  order,  dramage 
area  class,  or  Rosgen  Level  I stream  type 
designation. 

2.  The  differences  observed  in  the  taxa 
richness  and  modified  percent  mayfly- 
index  scores  of  Rosgen  Level  II  stream 
types  are  more  readily  attnbuted  to 
ecoregion.  rather  than  stream  reach 
channel  morphology  conditions 

3.  The  relationships  observed  between  both 
modified  EPT  and  percent  dominant  taxa 
index  scores  and  modified  Rosgen  stream 
type  appear  to  be  independent  of 
ecoregion  influences.  This  conclusion  is 
supported  by  the  fact  that  the  various 
modified  Rosgen  stream  types  were 
relatively  well  distnbuted  throughout  the 
different  ecoregions  included  in  the  study  . 
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Table  12.  Summarv  of  Statistical  Analyses  (ANOVA  and  t-test)  Between  Stream  Classification  Cariahles  and  Metric  Scores  of  Data  Set  3 
Sample  Sites  (Significant  differences  are  reported  at  the  95%  confidence  level  (p<().05).) 
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PERCENT  DOMINANT  TAXA 


ECOREGION 

(a) 


(b) 


Figure  13.  Box-and-lVhisker  Plots  of  (a)  Percent  Dominant  Taxa  and  (b)  Modified  Percent  Mayfly 
Index  Scores  of  Data  Set  3 Sample  Sites  by  Ecoregion 
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Figure  14.  Box-and-  Whisker  Plots  of  (a)  Taxa  Richness  of  Data  Set  3 Sample  Sites  and  (b)  Modified 
Percent  Mayfly  Index  Scores  of  Group  E (see  Table  12)  Sample  Sites  by  Rosgen  Level  II 
Stream  Type 
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(a) 


(b) 


Figure  15.  Box-and-Whisker  Plots  of  (a)  Modified  EFT  Index  Scores  of  Data  Set  3 Sample  Sites  and 
(h)  Percent  Dominant  Taxa  Index  Scores  of  Group  B (see  Table  12)  Sample  Sites  by 
Modified  Rosgen  Stream  Type 
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Figure  1 6.  Box-and-  Wlnsker  Plots  of  Modified  Hilsenhoff  Biotic  Index  Scores  of  Data  Set  3 Sample 
Sites  by  Fb'edominant  Channel  Material 
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For  example,  the  ANOVAs  that  disclosed 
significant  differences  in  the  modified 
EPT  and  percent  dominant  taxa  index 
scores  of  modified  Rosgen  stream  t\pes 
included  eight  sample  sites  from 
Ecoregion  62  (four  B3  streams,  one  C3 
stream,  and  three  C4  streams)  and  four 
sites  from  Ecoregion  60  (one  B3  stream, 
one  C3  stream,  and  two  C4  streams). 

4.  In  addition  to  being  an  effective  technique 
for  partitioning  variation  among  stream 
benthic  macro  invertebrate  communities, 
the  modified  Rosgen  technique  classifies 
streams  at  a level  of  detail  intermediate  to 
Rosgen  Levels  I and  II  and  has  much 
potential  to  be  used  as  a practical  cost- 
effective  method  for  classif-ing  streams. 


MANAGEMENT  IMPLICATIONS 
Stream  Classification 

The  data  included  in  this  study  were  useful  for 
identifying  several  relationships  between  stream 
classification  variables  and  the  benthic 

macromvertebrate  communities  of  HQ-CWFs.  At 
the  segment  spatial  scale  (Table  1).  significant 
differences  exist  in  both  the  percent  dominant  taxa 
and  modified  percent  ma\fly  index  scores  of 
sample  sites  based  on  ecoregion  designation.  At 
the  stream  reach  spatial  scale,  stream  order, 
drainage  area,  surface  water  slope,  and  Rosgen 
Level  I classification  are  not  effective 

classification  vanables.  However,  modified 

Rosgen  stream  t>pe  designation  is  effective  for 
partitioning  variation  among  stream  benthic 
macroinvertebrate  communities  of  HQ-CWTs  at 
the  stream  reach  spatial  scale.  Thus,  it  is 
concluded  that  the  combination  of  ecoregion  and 
modified  Rosgen  stream  t\pe  designation  is  the 
most  effective  stream  classification  scenano 
assessed  in  this  study. 

The  modified  Rosgen  technique  classifies 

streams  at  a level  of  detail  intermediate  to  Rosgen 
Levels  I and  II  and  has  much  potential  to  be  used 
as  a practical,  cost-effective  method  for 
classify  ing  streams.  Modified  Rosgen  stream  t\pe 


designation  is  generated  from  Rosgen  Level  I 
designation  enhanced  with  channel  matenal 
information.  Rosgen  (1996)  descnbes  Level  I 
classification  as  being  designed  to  be  a coarse- 
scale  geomorphic  charactenzation  that  can  be 
performed  using  information  taken  from 
topographic  and  landform  maps  and  aenal 
photographs,  followed  by  subsequent  venfication 
with  field  data  from  selected  reference  reaches. 
.Modified  Rosgen  stream  t>pe  designation  is  based 
pnmanly  on  the  stream  channel  morphology 
parameters  entrenchment  ratio,  width/depth  ratio, 
and  dominant  channel  materials.  Field  data 
pertaining  to  these  parameters  can  be  collected 
with  a minimum  of  effort,  and  it  is  quite  probable 
that,  with  expenence  and  field  calibration  and 
verification,  trained  personnel  could  sufficientlv 
estimate  these  parameters  in  the  field. 
Furthermore,  modified  Rosgen  stream  t\pe 
categories  could  be  broken  down  further  into 
more  detailed  stream  t\pe  classes  based  on 
additional  data  generated  m the  field  and/or  office. 
Examples  of  office-generated  data  mclude 
estimates  of  sinuosity  and  slope  taken  from  aenal 
photographs,  topographic  maps,  and/or  GIS. 

Areas  for  Additional  Research  and 
Development 

Developing  GIS  applications  for  accuratelv 
generating  and/or  managing  Rosgen  classification 
data  would  provide  water  resource  managers  with 
an  effective  tool  for  classify  ing  streams,  screening 
stream  reach  morphological  conditions,  organ- 
izing watershed  data,  developing  watershed 
management  guidelines,  and  effectively  targeting 
stream  and  watershed  restoration  activities.  Of 
particular  interest,  is  the  potential  use  of  Rosgen 
stream  classification  data  to  guide  watershed 
management  activities.  M\ers  and  Swanson 
(1992)  showed  that  various  Rosgen  stream  tspes 
responded  differently  to  riparian  grazing.  Rosuen 
(1994  and  1996)  expanded  on  the  findings  of 
Meyers  and  Swanson  (1992)  and  developed  a 
framework  for  characterizing  streams  with  respect 
to  their  sensitivity  to  disturbance.  reco\er\' 
potential,  sediment  supply,  streambank  erosion 
potential,  and  the  influence  of  vegetation  at 
controlling  stabiliw.  based  on  Rosgen  stream  tvpe 
designation.  This  mformation  can  be  used  to 


guide  land  and  watershed-specific  management 
activities,  based  on  the  Upies  of  streams  that  dram 
the  watershed.  Other  potential  uses  of  Rosgen 
classification  data  include  applications  related  to 
hydraulic  geometiy,  flow  resistance,  grazing 
management  techniques,  fish  habitat  improvement 
structure  suitability,  and  stream  channel 
restoration  design. 

Problems  Encountered 

Initially,  Rosgen  stream  t\pe  designations  of 
sample  sites  w'ere  made  using  the  channel  cross- 
sectional  dimensions  associated  with  bankfull 
stage  elevations  delineated  in  the  field.  These 
delineations  were  based  on  field  indicators  of 
bankfull  stage  such  as  floodplain  elevations,  the 
tops  of  depositional  features  such  as  point  bars, 
breaks  in  streambank  slope,  changes  m 
streambank  particle  size,  and  staining  of  rocks. 
However,  after  receiving  comments  pertaining  to 
some  of  the  bankfull  channel  dimensions  and 
cross  sections  included  m an  earlier  version  of  this 
publication,  the  bankfull  channel  dimensions  of 
each  sample  site  were  reviewed  and  compared  to 
curves  of  average  values  of  bankfull  channel 
dimensions  and  discharge  as  functions  of  drainage 
area  reported  in  Dunne  and  Leopold  (1978).  The 
outcome  of  this  review  suggested  that,  at  many  of 
the  sample  sites,  bankfull  stage  elevation  was 
substantially  underestimated  during  the  collection 
of  field  data.  Thus,  the  bankfull  stage  elevations 
of  all  sample  sites  were  reevaluated  using  the 
procedure  described  in  the  methods  section  of  this 
report  and  in  several  cases  were  raised  to  more 
realistic  elevations  and  subsequently  reclassified 
using  Rosgen’s  classification  procedure. 

Bankfull  stage  elevation  is  a critical 
component  of  the  Rosgen  classification  system 
and  unfortunately  it  also  is  quite  difficult  to 
accurately  delineate  in  the  field  (Johnson  and 
Heil,  1996).  At  the  onset  of  this  project,  the 
delineation  of  bankfull  stage  elevation  appeared  to 
be  a relatively  straightforward  task  that  could  be 
accomplished  in  the  field.  However,  this  did  not 
turn  out  to  be  the  case,  and  the  following 
recommendations  are  offered:  (1)  investigators 
should  perform  estimates  of  bankfull-stage 
elevation  at  appropriate  gaged  stream  reaches 


before  collecting  stream  reach  classification  data; 
and  (2)  investigators  should  obtain  or  develop 
regional  (or  ecoregional)  curves  of  average  values 
of  bankfull  channel  dimensions  and  discharge  as 
functions  of  drainage  area  before  using  Rosgen’s 
stream  classification  system.  These 

recommendations  will  improve  the  accuracy  of 
Rosgen  classification  designations  and  provide 
water  resource  managers  with  the  opportunity  to 
utilize  the  Rosgen  stream  reach  classification 
system  to  its  fullest  potential. 
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Appendix  A 

Channel  Cross-Sectional  and  Pebble-Count  Data 
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Figure  A 1.  Cross-Sectional  and  Pebble-Count  Data  From  Asaph  Run 
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Figure  A2.  Cross-Sectional  and  Pebble-Count  Data  From  Baldwin  Branch 
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Figure  A3.  Cross-Sectional  and  Pebble-Count  Data  From  Cork  Hill  Creek 
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Figure  A4.  Cross-Sectional  and  Pebble-Count  Data  From  Deep  Hollow 
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Figure  A5.  Cross-Sectional  and  Pebble-Count  Data  From  Dead  Woman  Hollow 
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Figure  A6.  Cross-Sectional  and  Pebble-Count  Data  From  East  Branch  Codorus  Creek 
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Figure  A 7.  Cross-Sectional  and  Pebble-Count  Data  From  East  Branch  Lackawanna  River 


47 


NUMBER  OF  PARTICLES 


O 

1- 

< 

> 

LU 

_l 

LU 


a: 

LU 


LL 


LU 

o 

cn 

LU 

Q. 

LU 

> 

1- 


13 

o 


CROSS-SECTIONAL  DATA 


109.5  - 


107.5  -I 


105.5  ^ 


103.5 


101.5 


99.5 


97.5  - 


95.5 


29 


Bankfull  Channel  Data 

Rosgen  Classification  Data 

Area  (ft^) 

44.0 

Entrenchment 

> 1.5 

Perimeter  (ft) 

29.3 

Width/Depth  Ratio 

17.0 

Width  (ft) 

27.2 

Sinuosity 

1.0 

Average  Depth  (ft) 

1.6 

Slope 

0.019 

Estimated  Discharge  (cfs) 

296.38 

Channel  Material 

Gravel 

Stream  Type 

B3c 

I 

49 

1 

69 

i 

89 

DISTANCE  (ft) 


PEBBLE-COUNT  DATA 


BED  SURFACE  MATERIAL  - PARTICLE  SIZE  CLASSES  (mm) 


Figure  A8.  Cross-Sectional  and  Pebble-Count  Data  From  East  Branch  Roaring  Brook 
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Figure  A9. 


Cross-Sectional  and  Pebble-Count  Data  From  Gazzam  Run 
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Figure  A 10,  Cross-Sectional  and  Pebble-Count  Data  From  Greenlee  Run 
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Figure  All.  Cross-Sectional  and  Pebble-Count  Data  From  Green  Run  (Lower) 


51 


NUMBER  OF  PARTICLES 


CUMULATIVE  PERCENT  FINER  ELEVATION  (ft) 


CROSS-SECTIONAL  DATA 


PEBBLE-COUNT  DATA 


BED  SURFACE  MATERIAL  - PARTICLE  SIZE  CLASSES  (mm) 


Figure  All.  Cross-Sectional  and  Pebble-Count  Data  From  Green  Run  (Upper) 
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Figure  A13.  Cross-Sectional  and  Pebble-Count  Data  From  Hammer  Creek  (Lower) 
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Figure  A14.  Cross-Sectional  and  Pebble-Count  Data  From  Hammer  Creek  (Upper) 
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Figure  A1 5.  Cross-Sectional  and  Pebble-Count  Data  From  Harveys  Creek 
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Figure  A 16.  Cross-Sectional  and  Pebble-Count  Data  From  Horse  Hollow 
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Figure  A 17.  Cross-Sectional  and  Pebble-Count  Data  From  Hutton  Run 
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Figure  A 18.  Cross-Sectional  and  Pebble-Count  Data  From  Iron  Run 
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Figure  A 19.  Cross-Sectional  and  Pebble-Count  Data  From  Lebo  Run 
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Figure  A20.  Cross-Sectional  and  Pebble-Count  Data  From  Little  Moores  Run 
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Figure  A2L  Cross-Sectional  and  Pebble-Count  Data  From  Letort  Spring  Run 
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Figure  A22.  Cross-Sectional  and  Pebble-Count  Data  From  Mountain  Creek 
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Figure  A23.  Cross-Sectional  and  Pebble-Count  Data  From  Murray  Run 
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Figure  A24.  Cross-Sectional  and  Pebble-Count  Data  From  Pequea  Creek 
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Figure  A25.  Cross-Sectional  and  Pebble-Count  Data  From  Rice  Run 
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Figure  A26.  Cross-Sectional  and  Pebble-Count  Data  From  Sixmile  Run 
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Figure  A27.  Cross-Sectional  and  Pebble-Count  Data  From  Standing  Stone  Creek 
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Figure  A28.  Cross-Sectional  and  Pebble-Count  Data  From  Tipton  Run 
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Figure  A29.  Cross-Sectional  and  Pebble-Count  Data  From  Toms  Run  (Lower) 
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Figure  A30.  Cross-Sectional  and  Pebble-Count  Data  From  Toms  Run  (Upper) 
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Figure  A3 1.  Cross-Sectional  and  Pebble-Count  Data  From  Unnamed  Tributary  to  East  Branch 
Codorus  Creek 
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Figure  A3 2.  Cross-Sectional  and  Pebble-Count  Data  From  Unnamed  Tributary  to  Harveys  Creek 
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Figure  A3 3.  Cross-Sectional  and  Pebble-Count  Data  From  Unnamed  Tributary  to  Roaring  Brook 


73 


NUMBER  OF  PARTICLES 


CUMULATIVE  PERCENT  FINER  ELEVATION  (ft) 


CROSS-SECTIONAL  DATA 


PEBBLE-COUNT  DATA 


^ 04 

BED  SURFACE  MATERIAL  - PARTICLE  SIZE  CLASSES  (mm) 


Figure  A34.  Cross-Sectional  and  Pebble-Count  Data  From  Willow  Run 
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Fisure  A35.  Cross-Sectional  and  Pebble-Count  Data  From  White  Deer  Creek 
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Class:  Order 

Family 

Family/Genus 

Tolerance 

Trophic  Guild 

ASPH 

BALD 

CORK  1 DEEP 

Insecta:  Coleoptera 

ilmidae 

Optiosennis 

4 

SC 

-) 

Oulimnms 

5 

SC 

1 

Promoresia 

2 

SC 

Stenelmis 

5 

SC 

1 

Psephenidae 

Psephemis 

4 

SC 

1 

Ectoprm 

5 

SC 

1 

1 

Ptilodactylidae 

Anchytarsus 

5 

SH 

Diptera 

\thericidae 

Athenx 

2 

P 

5 

Ceratoposonidae 

Bezzia 

6 

P 

1 

Chironomidae 

Chironomidae 

7 

CG 

2 

14 

“1 

1 

Empididac 

Hetnerodromia 

6 

P 

Simuliidae 

Simuliidae 

6 

FC 

Tabanidae 

Tabamis 

5 

P 

1 

Tipulidae 

Antocha 

3 

CG 

2 

Dicranota 

3 

P 

1 

1 

H exatom  a 

2 

P 

3 

2 

2 

Lunoma 

6 

SH 

Tipnla 

4 

SH 

Ephemeroptera 

Ameletidae 

Aineletus 

0 

CG 

Baetidae 

Acentrella 

4 

CG 

1 

Baetis 

6 

CG 

12 

2 

Ceniroptilum 

2 

CG 

Cloeon 

4 

CG 

Baetiscidae 

Baetisca 

4 

CG 

Ephemerellidae 

Ephemerella 

1 

SC 

27 

6 

Eiirylophella 

4 

SC 

Serratella 

2 

SC 

7 

Ephemeridae 

Ephemera 

2 

CG 

Heptageniidae 

Epeoms 

0 

CG 

Heplagema 

4 

SC 

2 

7 

Leucrocuta 

1 

SC 

Rhithro^ena 

0 

CG 

15 

Stenacron 

4 

SC 

Stenonema 

3 

SC 

7 

8 

3 

2 

Isonychiidae 

honvchia 

2 

FC 

Leptophlebiidae 

Paraleptophlebia 

1 

CG 

4 

1 

6 

1 

Tricorythidae 

Tricorythodes 

4 

CG 

Megaloptera 

Corydalidae 

Ni^roma 

2 

P 

2 

2 

Sialidae 

Sialis 

4 

P 

Odonata 

Aeshnidae 

Boyena 

2 

P 

1 

Gomphidae 

Lanthus 

5 

P 

1 

1 

Stylogomphus  albistylus 

4 

P 

Plecoptera 

Chloroperlidae 

Alloperla 

0 

CG 

3 

10 

2 

Haploperla 

0 

P 

2 

Sweltsa 

0 

P 

2 

Leuctridae 

Leuctra 

0 

SH 

Peltoperiidae 

Peltoperla 

2 

SH 

Tatlaperla 

0 

SH 

1 

4 

4 

Perlidae 

Acroneuna 

0 

P 

2 

1 

Asnetina 

2 

P 

1 

Eccoptura 

2 

P 

3 

Para^netina 

1 

P 

2 

Perlodidae 

Cultus 

2 

P 

Diploperla 

2 

P 

Hydroperla 

1 

P 

5 

8 

2 

Isoperla 

2 

P 

1 

4 

2 

Rememis 

2 

P 

Pteronarc\idae 

Pteronarcys 

0 

SH 

4 

16 

Taeniopterysidae 

Strophopteryx 

3 

SH 

Taemopteryx 

2 

SH 

4 
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Class:  Order 

Family 

Family/Genus 

Tolerance 

Trophic  Guild 

ASPH 

BALD 

GORKI  DEEP 

Trichoptera 

Brachyccntridae 

Brachyceninis 

1 

FC 

Glossosomatidae 

Asapetus 

0 

SC 

Glossosoma 

0 

sc 

1 

1 

Protoplila 

1 

sc 

Hydropsychidae 

CheiimatoDsyche 

5 

FC 

5 

2 

11 

3 

Dipleclrona 

0 

FC 

6 

Homoplectra 

0 

FC 

Hydropsyche 

4 

FC 

7 

6 

6 

12 

Potamyta  flava 

5 

FC 

1 

3 

Lepidostomatidae 

Theliopsyche 

1 

SH 

Odontoceridae 

Psilotreta 

0 

SC 

Philopotamidae 

Chimarra 

4 

FC 

Dolophilodes 

0 

FC 

21 

17 

21 

24 

Polycentropodidae 

Neureclipsis 

7 

FC 

2 

Polycentropus 

6 

FC 

Psychomyiidae 

Lvpe  diversa 

2 

CG 

Rhyacophilidae 

Rhyacop’mla 

1 

P 

3 

4 

Rh\Tichobdellida 

Glossiphonidae 

Batrachobdella 

6 

P 

Crustacea;  Amphipoda 

Gamma  ridae 

Gammanis 

6 

SH 

Decapoda 

Cambaridae 

Cambarus 

6 

CG 

1 

Isopoda 

Asellidae 

Caecidotea 

8 

SH 

Lircetis 

8 

CG 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 

6 

SC 
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Class:  Order 

Family 

Family/Genus  | 

DWMN  1 

EBCC  1 

EBLR 

EBRB 

GAZZ  1 

GRLE 

Insecta:  Coleoptera 

Llmidae 

Optioserviis 

6 

8 1 

4 

1 

13  1 

Ouhmnnis 

^ i 

Promoresia 

1 

2 1 

Stenelmis 

’sephenidae 

Psevhenus 

1 

2 

4 

Ectopna 

1 

1 

Ptilod  jK>'lidae 

4nchyrarsus 

Diptera 

\thericidae 

-Ithenx 

4 

2 

-1 

Ceratopogonidae 

Be::ia 

2 

Chironomidae 

Chironomidae 

1 

5 

4 

4 

14 

8 

Empididae 

Hemerodromia 

Simuliidae 

Simuliidae 

1 

Tabanidae 

Tabanus 

Fipulidae 

Antocha 

Dicranota 

1 

1 

1 

Hexaloma 

5 

6 

3 

1 

Limonia 

Tipula 

Ephemeroptera 

.Ameletidae 

Ameletus 

Baetidae 

Acentrella 

8 

2 

Baetis 

5 

Centropiilum 

3 

Cloeon 

Baetiscidae 

Baelisca 

Ephemerellidae 

Ephemerella 

19 

1 

Eurylophella 

Serratella 

3 

Ephemeridae 

Ephemera 

Heptageniidae 

Epeonis 

Hepta^ema 

Leucrocuta 

Rhithro°ena 

Stenacron 

Stenonema 

6 

15 

15 

3 

8 

Isonychiidae 

Isonychia 

8 

6 

Leptophlebiidae 

Paraleptophlebia 

8 

5 

6 

Tricorythidae 

Tncorythodes 

Megaloptera 

Corydalidae 

Nigronia 

3 

1 

6 

2 

2 

Sialidae 

Sialis 

1 

Odonata 

Aeshnidae 

Bovena 

Gomphidae 

Lanthus 

1 

2 

1 

Srvlooomphus  aibisrylus 

1 

Plecoptera 

Chloroperiidae 

Alloperla 

9 

Haploperla 

Swehsa 

Leuctridae 

Leuctra 

Peltoperlidae 

Peltoperla 

Tallaperla 

4 

Perlidae 

Acroneuria 

7 

1 

6 

3 

8 

A^netina 

5 

Eccopmra 

4 

1 

ParaoneUna 

Perlodidae 

Cultiis 

Diploperla 

Hydroperla 

Isoperla 

1 

1 

1 

Rememis 

Pteronarcyidae 

Pteronarcvs 

1 

1 

T aeniopterygidae 

Slrophopterx'x 

1 TaenioDteryx 

6 

10 

1 
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Class:  Order  | Family 

Famiiy/Genus 

DWMN 

EBCC 

EBLR 

EBRB 

GAZZ 

GRLE 

Trichoptera 

Brachycentridae 

Brachycenlnis 

Glossosomatidae 

Asapetus 

1 

11 

Glossosoma 

1 

Protopula 

3 

Hydropsychidae 

Cheumatopsyche 

4 

13 

2 

1 

7 

2 

Diplectrona 

3 

8 

2 

Homoplectra 

Hydropsyche 

30 

21 

8 

9 

22 

Potamyia  flava 

2 

25 

Lepidostomatidae 

Theliopsyche 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

1 

Dolophilodes 

2 

2 

7 

6 

4 

25 

Poiycentropodidae 

Neureclipsts 

3 

1 

Polycentropus 

Psychomyiidae 

Lvpe  diver sa 

Rhyacophiiidae 

Rhvacophila 

1 

Rhynchobdellida 

Glossiphonidae 

Batrachob  della 

Crustacea:  Amphipoda 

Ganunaridae 

Gammanis 

Decapoda 

Cambaridae 

Cambanis 

2 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Levtoxis 
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Class:  Order  | Family 

Famiiy/Genus 

GRNL  1 GRNU  | HAML 

HAMU  1 HARV  | HORS 

Insecta;  Coleoptera  1 

-Imidae 

Dpuosennis 

1 

3 

22 

8 

6 

Oulimmus 

^romoresia 

1 

Stenebms 

4 

3 

I 

’sephenidae 

Psephenus 

J 

Ectopna 

I 

’tilodact>'lidae 

Anchytarsus 

Diptera 

Vthericidae 

Atherix 

( 

reratoposonidae 

Bezzia 

( 

ITiironomidae 

Chironomidae 

7 

6 

3 

34 

11 

^mpididae 

Hemerodrotnia 

Miniiliidae 

Simuliidae 

1 

2 

2 

Tabanidae 

Tabamis 

1 

ripulidae 

Antocha 

4 

Dicranota 

Hexatoma 

1 

1 

2 

Limoma 

Tipula 

1 

Ephemeroptera 

\meleridae 

Ameletus 

1 

Jaetidae 

Acentrella 

3 

Baens 

12 

8 

15 

Centropuhim 

Cloeon 

Baetiscidae 

Baetisca 

Ephemerellidae 

Ephemerella 

Eurylophella 

2 

Serratella 

3 

1 

Ephetneridae 

Ephemera 

1 

Heptageniidae 

Epeonis 

Hepta^ema 

Leucrocuta 

Rhithro^ena 

Stenacron 

Stenonema 

4 

13 

10 

21 

Isonychiidae 

Isonychia 

17 

Leptophlebiidae 

Paraleptophlebia 

8 

6 

1 

T ricorythidae 

Tncorythodes 

1 

MegaJoptera 

Corydalidae 

Nigronia 

1 

2 

1 

9 

Sialidae 

Siahs 

Odonata 

Aeshnidae 

Boyena 

1 

2 

Gomphidae 

Lanthiis 

1 

Stylozomphus  albistylus 

1 

Plecoptera 

Chloroperlidae 

Alloperla 

2 

Haploperla 

Sweltsa 

Lcuctridae 

Leuctra 

Peltoperlidae 

Pelloperla 

Tallaperla 

1 

Perlidae 

Acroneuna 

6 

Aznetina 

Eccoptura 

2 

Paraznetina 

10 

Perlodidae 

Cultiis 

Diploverla 

3 

Hydroperia 

Isoperla 

*) 

Rememts 

Pteronarcyidae 

Pteronarcys 

Taeniopterysidae 

Strouhopieryx 

Taemopleryx 

3 
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Class;  Order 

Family 

Family/Genus  | GRNL 

GRNU 

HAML 

HAMU 

HARV 

HORS 

Trichoptera 

Brachycentridae 

Brachycentrus 

Glossosomatidae 

Agapetus 

Glossosoma 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

17 

15 

9 

Diplectrona 

4 

4 

1 

Homoplectra 

Hydropsyche 

6 

3 

45 

23 

6 

Potamyia  flava 

8 

16 

Lepidostomatidac 

Theliopsyche 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chunarra 

1 

21 

Dolophilodes 

15 

11 

Polycentropodidae 

Neureclipsis 

Polycentropus 

Psychomyiidae 

Lvpe  diversa 

2 

Rhyacophilidae 

Rhyacophda 

1 

4 

Rliynchobdellida 

Glossiphonidae 

Batrachob  della 

Crustacea:  .Vmphipoda 

Gammaridae 

Gammarus 

18 

1 

Decapoda 

Cambaridae 

Cambarus 

1 

2 

Isopoda 

Asellidae 

Caecidotea 

2 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order 

Family 

Family/Genus  | 

HUTT 

IRON 

LEBO 

LMOR 

LTRT  1 

MNTN 

Insecta:  Coleoptem 

Llmidae 

Optwservus 

2 

1 

6 

Dulimnms 

1 

3 

1 

Promoresia 

4 

Stemhms 

1 

’sephenidae 

Psephenus 

1) 

Ectopna 

2 

1 

?tilodact>lidae 

inchvtarsus 

Diptera 

\thericidae 

■Ithenx 

20 

1 

"eratopogonidae 

Bezzia 

Chironomidae 

Chironomidae 

1 

3 

2 

7 

15 

Empididae 

Hetnerodromia 

Sunuliidae 

Simuliidae 

1 

1 

Tabanidae 

Tabanus 

Tipulidae 

-intocha 

Dicranota 

2 

1 

2 

Hexatoma 

9 

5 

10 

Limonia 

Tipida 

3 

1 

1 

1 

1 

Ephemeroptera 

Ameletidae 

Ameletus 

Baeddae 

Acentrella 

Baetis 

26 

8 

2 

1 

Centroptilum 

Cloeon 

Baetiscidae 

Baetisca 

Ephemerellidae 

Epheinerella 

15 

1 

1 

Eurylophella 

Serratella 

4 

Ephemeridae 

Ephemera 

Heptageniidae 

Epeortis 

3 

Heptas,ema 

3 

13 

Leucrocuta 

4 

Rhithro^ena 

Stenacron 

Stenonema 

3 

1 

4 

8 

Isonychiidae 

honychia 

Leptophlebiidae 

Paraleptophlebia 

1 

Tricorythidae 

Tncorythodes 

Megaloptera 

Corydalidae 

Ni^roma 

2 

1 

Sialidae 

Slabs 

Odonata 

Aeshnidae 

Boyeria 

Gomphidae 

Lanthus 

5 

Stylo^omphus  albistylus 

Plecoptera 

Chloroperlidae 

Alloperla 

Haploperla 

Sweltsa 

Leuctridae 

Leuctra 

Peltoperiidae 

Peltoperla 

Tallaperla 

3 

3 

8 

Perlidae 

Acroneuna 

A^netina 

2 

Eccoptura 

6 

1 

Paraonetma 

Perlodidae 

Cultus 

Diploperla 

Hvdroperla 

2 

4 

1 

Isoperla 

6 

1 

3 

4 

Remenus 

Pteronarcyidae 

Pieronarcys 

3 

4 

1 

T aeniopterygidae 

Strophopteryx 

1 

TaemoDteryx 
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Class;  Order 

Family 

Famiiy/Genus 

HUTT 

IRON 

LEBO 

LMOR 

LTRT 

MNTN 

Trichoptera 

Brachycentridae 

Brachycentnis 

1 

1 

Glossosomatidae 

Agapetus 

Glossosoma 

2 

1 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

17 

Diplectrona 

11 

1 

3 

Homoplectra 

Hydropsyche 

15 

4 

9 

10 

Potamyia  flava 

15 

6 

8 

Lepidostomatidae 

Theliopsvche 

6 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

Dolophilodes 

10 

13 

5 

2 

Polycentropodidae 

Neureclipsis 

Polycentropiis 

Psychomyiidae 

Lype  diversa 

Rhyacophilidae 

Rhyacophila 

2 

3 

1 

Rhynchobdellida 

Giossiphonidae 

Batrachob  della 

Crustacea:  Amphipoda 

Gammaridae 

Gammarus 

63 

Decapoda 

Cambaridae 

Cambarus 

1 

1 

1 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

37 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order  | Family  | Family/Genus  | MURY 

PEQU 

RICE  1 

SIXM  1 

STSN  1 TIPT 

Insecta:  Coleoptera 

Slmidae 

Dvttoser\nis 

1 

3 

1 

Dulimmiis 

1 

i 

1 

°romoresia 

2 

Stenelmis 

’sephenidae 

Psephenns 

1 

1 

J 

Ectopna 

2 

3 

’tilodactA’lidae 

Anchvtarsus 

Diptera 

Xthericidae 

Atherix 

leratopogonidae 

Bezzia 

Chironomidae 

Chironomidae 

12 

1 

3 

5 

6 

8 

Smpididae 

Hemerodromia 

1 

Simuliidae 

Simuliidae 

13 

4 

1 

Tabanidae 

Tabanus 

Tipulidae 

Antocha 

6 

Dicranota 

Hexatoma 

2 

6 

2 

1 

Limonia 

Tipula 

Ephemeroptera 

.\meletidae 

Ameletus 

Baetidae 

Acentrella 

1 

2 

Baetis 

4 

2 

21 

Centroptilum 

Cloeon 

4 

Baetiscidae 

Baetisca 

1 

Ephemerellidae 

Ephemerella 

2 

Eurylophella 

Serratella 

1 

Ephemeridae 

Ephemera 

Heptageniidae 

Epeorus 

Hepta^ema 

Leucrocuta 

4 

Rhithrogena 

Stenacron 

Stenonema 

13 

4 

8 

2 

3 

Isonychiidae 

Isonychia 

1 

3 

Leptophlebiidae 

Paraleptophlebia 

2 

1 

1 

Tricorythidae 

Tncorythodes 

Megaloptera 

Corydalidae 

Ni^roma 

1 

5 

3 

1 

Sialidae 

Stalls 

1 

Odonata 

Aeshnidae 

Boyeria 

Gomphidae 

Lanthus 

4 

1 

Stylooomphus  albistyhis 

1 

1 

Plecoptera 

Chloroperlidae 

Alloperla 

4 

4 

Haploperla 

Sweltsa 

Leuctridae 

Leuctra 

4 

Peltoperlidae 

Peltoperla 

2 

Tallaperla 

1 

1 

Perlidae 

Acroneuna 

7 

8 

A°netina 

Eccoptura 

9 

1 

Paragnetina 

2 

Perlodidae 

Cultus 

Diploperla 

Hvdroperla 

Isoperla 

9 

5 

Rememis 

Pteronarcyidae 

Pteronarcys 

19 

3 

1 

Taeniopterygidae 

Strophopteryx 

Taemopteryx 
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Class:  Order  | Family 

Family/Genus  1 MURY 

PEQU 

RICE 

SIXM 

STSIM 

TIPT 

Trichoptera 

Brachycentridae 

Brachycentnis 

Glossosomatidae 

A^apetus 

Glossosoma 

1 

Protopttla 

Hydropsychidae 

Cheumatopsyche 

1 

2 

4 

1 

14 

Diplectrona 

3 

1 

9 

Homoplectra 

Hydropsyche 

46 

19 

4 

6 

7 

Potamyia  flava 

3 

7 

3 

12 

3 

Lepidostomatidae 

Theliopsyche 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

6 

Dolophilodes 

21 

4 

5 

21 

24 

41 

Polycentropodidae 

Neiireclipsis 

Polyceniropiis 

2 

Psychomyiidae 

Lype  diversa 

2 

Rhyacophilidae 

Rhyacophila 

3 

4 

10 

10 

3 

RhvTichobdellida 

Glossiphonidae 

Balrachobdella 

1 

Crustacea:  Amphipoda 

Gammaridae 

Gammarus 

Decapoda 

Cambaridae 

Cambanis 

1 

1 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order 

Family  1 Family/Genus 

TOML  1 TOMU  1 UEBC 

UHAR  1 

URRB  1 WILL 

WTDR 

Insecta:  Coleoptera 

ilmidae 

Jutioserviis 

4 

17 

1 

10 

Oulimmus 

2 

3 

J 

Promoresia 

4 

1 

10 

Stenelmis 

’sephenidae 

Psephenus 

1 

1 

J 

Ectovnc 

2 

^tilodact>iidae 

■Inchytarsus 

1 

1 

1 

Diptera 

A.thericidae 

Atherix 

10 

reratopogonidae 

Bezzia 

rhironomidae 

Chironomidae 

4 

4 

15 

6 

9 

Empididae 

Hemerodromia 

Simuliidae 

Simuliidae 

3 

18 

Tabanidae 

Tabamis 

Tipulidae 

Antocha 

1 

Dicranota 

6 

1 

1 

8 

Hexatoma 

2 

3 

3 

3 

Lunoma 

1 

Tivula 

6 

1 

Ephemeroptera 

Ameletidae 

Ameletiis 

Baetidae 

AcenP-ella 

Baetis 

2 

2 

Centropnlum 

2 

Cloeon 

2 

Baetiscidae 

Baensca 

Ephemerellidae 

Ephemerella 

7 

4 

Eun'lophella 

Serratella 

Ephemeridae 

Ephemera 

1 

Heptageniidae 

Epeorus 

1 

1 

Hepta^enia 

Leucrocuta 

Rhithrogera 

Stenacron 

1 

Stenonema 

1 

5 

19 

Isonychiidae 

honvchia 

30 

Leptophlebiidae 

Paraleptophlebta 

21 

5 

T ricorythidae 

Tncorythodes 

Megaloptera 

Corydalidae 

Nizronia 

7 

6 

2 

1 

5 

9 

Sialjdae 

Sialis 

Odonata 

Aeshnidae 

Bovena 

1 

Gomphidae 

Lanthus 

2 

St\iop,omphus  albistylus 

1 

Plecoptera 

Chloroperlidae 

Alloperla 

2 

Haploperla 

Sweltsa 

Leuctridae 

Leuclra 

2 

Peltoperlidae 

Peltoperla 

Tallaperla 

1 

3 

Perlidae 

Acrorjeuna 

1 

3 

6 

12 

1 

3 

Agnenna 

Eccoplura 

5 

3 

Paragnenna 

Perlodidae 

Culna 

1 

Dwloperla 

Hx’droperla 

6 

Isoperla 

1 

Remenus 

Pteronarcyidae 

Pteronarcys 

1 

1 

T aeniopterygidae 

Strophopteryx 

1 

7 

Taemopteryx 
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Class:  Order 

Family 

Family/Genus 

TOML 

TOMU 

UEBC 

UHAR 

URRB 

WILL 

WTDR 

Trichoptera 

Brachycentridae 

Brachycentnis 

4 

Glossosomatidae 

Agapetus 

Glossosoma 

4 

5 

1 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

4 

4 

7 

14 

7 

5 

9 

Diplectrona 

3 

1 

2 

2 

Homoplectra 

4 

Hydropsyche 

2 

23 

5 

3 

3 

4 

Potamyia  flava 

4 

11 

2 

Lepidostomatidae 

Theliopsyche 

1 

Odontoceridae 

Psilotreta 

3 

Philopotamidae 

Chiinarra 

19 

1 

4 

Dolophilodes 

2 

1 

13 

11 

16 

2 

16 

Polycentropodidae 

Neureclipsis 

Polycentropus 

Psychomyiidae 

type  diverse 

3 

Rhyacophilidae 

Rhyacophila 

1 

2 

1 

2 

5 

Rhynchobdellida 

Glossiphonidae 

Batrachobdella 

1 

Crustacea:  Amphipoda 

Ganunaridae 

Gammarus 

Decapoda 

Cambaridae 

Cambanis 

2 

1 

Isopoda 

Asellidae 

Caecidotea 

Urceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 

1 
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